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1. Introduction 

The B c meson provides a unique probe of both strong and weak interactions. Unlike 
quarkonium systems which can decay strongly and electromagnetically, the B c can only 
decay weakly and thus is relatively long-lived. The c and b quark lifetimes are similar 
because of the small mixing angle V c b, so B c decay proceeds through either quark at 
comparable rates. In this work, we study semileptonic weak decay of the B c meson by 
exploiting heavy quark spin symmetry. We will compute the weak decay amplitude for 
B c — > B s in terms of the B s meson decay constant fs s - The two are related because the 
B c can be treated as a pointlike meson in the limit that the b and c quark masses are much 
larger than Aqcd- The weak decay of the c quark produces a state which has a b and s 
quark at the same point in space (to lowest order in an operator product expansion); the 
amplitude for this state to turn into a B s is fs a • A similar argument allows us to compute 
the amplitudes for B c — > D°, D*° and B c — ► r] c , J/ifj semileptonic decay. 

The standard application of heavy quark symmetry is to hadrons containing a single 
heavy quark. We must take some care in defining the heavy quark effective theory when 
dealing with a system with two heavy quarks. It is well known that the static theory gives 
rise to severe infrared divergences for diagrams involving two heavy quarks with the same 
velocity [0. These divergences are regulated by the kinetic energy term /iQ(D 2 /2mg)/ig; 
even though this term is higher order in 1/mg, it may not be neglected in the uiq — > oo 
limit. The kinetic energy term is different for b and c quarks, and breaks the heavy flavor 
symmetry. Physically, this is just the statement that the dynamics of heavy-heavy bound 
states is determined by balancing the kinetic and potential energies of the quarks; the T is 
not the same size as the J/ip, and we cannot use the heavy flavor symmetry to relate these 
two states. The kinetic term in the effective Lagrangian breaks the heavy flavor symmetry, 
but it does not break the heavy quark spin symmetry. Thus we can still derive relations 
for hadrons with two heavy quarks using heavy quark spin symmetry. 

2. Spin Symmetry 

The invariance of the effective Lagrangian under individual spin rotations on the b and 
c quarks allows us to relate the form factors for vector and axial vector currents between 
the B c and various pseudoscalar and vector mesons in the same way as for heavy-light 
systems [0. 
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Let us first consider the semileptonic decays B c — > B s (B°)e + u or S c — > 5*(I?* )e + f in 
which the the B c decays into a -B meson containing a light quark. These decays correspond 
to the semileptonic weak decay of the charm quark into a light s or d quark. Since the mass 
of the b quark is much greater than that of the c quark, the energy released in the decay of 
the c quark is much smaller than m&, and the b quark is not deflected. Thus the velocity of 
the final meson is the same as the velocity of the initial meson. The initial momentum of 
the B c is p^ — mB c v^, and the final momentum of the B a is p'^ — tub^ + q^ 1 , where q is 
a small residual momentum. The final B a is on shell, so q ■ v = 0(1/ tub)- The momentum 
transfer to the lepton system is 

k» =p» -p'» = (m Bc -m Ba )v^ - q*. (2.1) 

Heavy quark spin symmetry implies that the pseudoscalar B c meson is degenerate with 
the vector B* meson. The consequences of spin symmetry for hadronic matrix elements 
may be derived using the commutation relations of Isgur and Wise , or more compactly 
using the well-known trace formalism M . The lowest- lying be bound states are represented 
by a 4 x 4 matrix 

H ^ = (A+£) [B ,^_ Be7s] (A^l , (2 . 2) 

where B c and B* annihilate pseudoscalar and vector meson bc bound states of velocity v, 
respectively. A subscript v on the heavy meson fields has been suppressed. Under spin 
symmetries on the heavy quark and antiquark, the heavy meson field transforms as 

An analogous definition for _£Z"( CC ) describes the (t] c , J/ip) spin multiplet. The spin multiplet 
for the B a and B* is given by 

H { ? = [Bl^, - B alb ] ilzJ& , (2.4) 

where the subscript a = 1, 2, 3 (or u, d, s) is an SU(3) V flavor index. The field H ( a b) is a 
doublet under heavy quark spin symmetry and a 3 under flavor SU(3)v symmetry [|J, 

H®^(lJH®) a Sl. (2.5) 

Note that the pseudoscalar and vector meson fields B and B* have dimension 3/2 because 
they contain factors of ^Jtob and y/rriB* relative to the standard normalization for scalar 
and vector fields. 
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The amplitudes for semileptonic B c decay to B a and B* are determined by the matrix 
elements of the weak hadronic current <? a 7^(l — 75 )c between the meson states. The most 
general form for the matrix element of the current which respects the heavy quark spin 
symmetry is 

Bi*\v,q\q a Tc\B c ,v) = -y/m Bc m Ba Tr (h ( ? n(v,a q) T H^j , (2.6) 
where 

Q,(v, a q) = Qi(a q) + a Q O 2 (a <?) 4, (2.7) 

is the most general Dirac matrix that can be written in terms of the vectors q and v (recall 
that q ■ v=0). Terms with factors of f can be omitted because of the identities 

i>H&> = H^ cl \ H^j = -H^ b \ jH a = H a , H a j> = -H a . (2.8) 

Note that the factor of V multiplying ij( cb ) in Eq. ( |2.6| ) is required by the heavy quark 
spin symmetry on the c quark. The radius of the B c meson, ao, is the typical scale for 
the variation of the form factors (as will be shown in the next section). For a Coulomb 
bound state, a^ 1 ~ a s (a^ -1 ) m c ; the linear confining term in the potential makes the state 
somewhat smaller than this estimate. In our case, ao(B c ) ~ a$(J j^) — (500 MeV) -1 |J. 
Explicit evaluation of Eq. (|2.6| ) gives 

(B a , v, q\Vft\B e , v) = v / 2m Bc 2m Ba [fii + a 2 qyj , 

(B*, v, q\V^\B c , v)=- i v /2m Bc 2m B * a 2 e^p e w q a v p , (2.9) 

(B*, v, q\A^\B cl v) = v /2m Bc 2m B * [fix e* + a 2 e* ■ q uj , 

where and refer to the vector and axial vector currents Q a 7 At c and ? a 7 At 75C, respec- 
tively, and e M is the polarization vector of the B*. The form factor Q 2 is irrelevant for 
semileptonic B c — > B S (B°) decay because the contribution of Q 2 to the decay amplitude 
will be proportional to the lepton mass. In addition, Q 2 does not contribute to decay 
amplitudes at zero recoil, q = 0. Note that the dimensionless functions Qi(aoq) are inde- 
pendent of the light quark flavor index in the SU(3)v limit. Thus, the ratio of KM mixing 
angles |Ks/Kd| can be extracted from comparison of B c semileptonic decay to B s , B* and 
B°, B*°. Leading ,S£7(3) ^-violating light quark flavor dependence of the form factors may 
be estimated in chiral perturbation theory 0. 

A similar analysis applies to the decays B c — > D° and B c — > D*° in which the b 
quark decays to a u. In this case, however, the light antiquark will typically recoil with 
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momentum comparable to or larger than the c quark mass. In order for the final meson 
to be bound, there must be a correspondingly large momentum transfer to the spectator 
c quark, and the effective theory breaks down. Higher dimensional operators which we 
have neglected will be of order q/m c and will dominate for large momentum transfer to 
the light degrees of freedom. Our results are thus valid only for q <C m c , i.e. near the zero 
recoil point. With this caveat, the analysis proceeds exactly as before. The amplitudes 
can be written in terms of two invariant functions £1(00(7) and £2(000), 

(-D , v, q\Vn\B c , v) = v /2m Bc 2ro D [£1 + a £ 2 o M ] , 

(D*°, v, q\V^\B Cl v)=- iy/2m B 2m D * o £ 2 e^ a pe w q a v p , (2.10) 
(D*°, v, q\A^\B c , v) = ^j2m B 2m D * [£ x e* + a £ 2 e* • q « M ] • 

Measurement of this decay provides a means of determining the KM angle \V u b\. 

Finally, we analyze the semileptonic decays of the B c to the charmonium mesons r\ c 
and ip. Once again, the momentum transfer to the produced c quark may be large and our 
results are only valid near the zero recoil point. In this case, however, there is an additional 
spin symmetry of the produced antiquark, which forbids a form factor proportional to ^. 
Thus, the matrix elements for the semileptonic decay of B c to r\ c and ip near zero recoil 
are determined by a single function A(ooq): 

(rjc, v, q\Vn\B c , v) = y /2m Bc 2m r]c A v M , 
(-0, v, q\A^\B c , v) =^2m Bc 2m^ A e* , 

where and refer to the vector and axial vector currents Try^c and 67^750, respectively, 
and e M is the polarization vector of the ip. 

3. The Scale of Variation of Form Factors 

The invariant tensors in Eq. ( |2.6| ) are multiplied by dimensionless form factors Qi(aoq). 
In this section, we explain in greater detail why the scale of variation of the form factors is 
set by ao, the radius of the B c bound state. Before discussing the case of interest, it is useful 
to first consider the scale of variation of form factors for two different circumstances — the 
matrix elements for a heavy quark current between two heavy-light mesons and the matrix 
elements for a light quark current between two heavy-light mesons. Let us first analyze 
the well-known example of the matrix elements of a heavy quark current between meson 
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states containing a single heavy quark. For concreteness, consider the matrix elements of 
the current bTb between B and B* mesons. The matrix elements can be evaluated using 
the trace formalism, 

(B, v'\bTb \B,v)=m B Tr (h^) Sl(v, v) T^j , (3.1) 

where v and v' are the velocities of the initial and final meson fields, respectively. The Dirac 
matrix coupling the heavy quark spin indices in the trace must be V by the spin symmetry. 
The matrix O coupling the light quark indices is not constrained, and is the most general 
possible Dirac matrix that can be constructed out of v and v' . From Eq. (|2.8| ), it follows 
that there is only one possible invariant for O, the Isgur-Wise function £(v ■ v'). This 
dimensionless nonperturbative function is the form factor for the light degrees of freedom 
in the heavy-light mesons. The light degrees of freedom of a B meson with definite velocity 
v carry a momentum which is typically of order Aqcd^- Thus, the momentum transfer 
to the light degrees of freedom in the above transition is of order Aqcd(i> — v'). Since 
hadronic form factors for the light degrees of freedom vary on the momentum scale Aqcd 5 
the variation of the Isgur-Wise function is controlled by Aqcd(^ — v')/Aqcd = v — v'. 
Thus the Isgur-Wise function £(v ■ v') varies on the scale over which v ■ v' changes by order 
one.0 

The scale of variation of form factors is different for the matrix elements of a light quark 
current between the same states. In the following, we will assume that the momentum 
transfer of the transition is small compared with the mass of the heavy quark, so that the 
velocity of the B meson is not changed. The matrix elements of the light quark current 
dq^d are given by 

(B, v, q\ dTd \B, v, 0) = m B Tr (h^ n r (v, q) H®j , (3.2) 

where the states are described by both a velocity v and a residual momentum q || . The 
Dirac matrix coupling the heavy quark indices is the identity matrix because of the heavy 
quark spin symmetry. The Dirac matrix Q-p coupling the light quark indices is the most 

1 The B* — > Bj electromagnetic transition amplitude due to the b quark current 67^6 can be 
computed from Eq. ( [3.1|) with the substitution T = 7^. The amplitude vanishes; thus the b quark 
magnetic moment transition is a l/m& effect. 
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general Dirac matrix which transforms as V under Lorentz transformations. For V = 7^, 
we fincH 

Or = O^^/Aqcd) + 2 (q7Aqcd) ct^/Aqcd, (3.3) 

where we have used current conservation in writing Eq. ( |3.3|) . In contrast to the first 
example, the momentum transfer to the light degrees of freedom is q, so that the scale 
of variation of the form factors is now q/Aqcd- The form factors in Eq. ( |3.3|) have a 
variation on the scale p ■ p' ~ Aq CD (or v ■ v' ~ Aq CD /m^) instead of the scale v ■ v' ~ 1 
for the Isgur-Wise function £(v ■ v'). 

We now consider the scale of variation for the form factors found in Sect. 2. The 
B c — > B s transition amplitude is an example of a matrix element of an operator containing 
both a light quark and a heavy quark. The matrix elements of the current sTc are given 
by 

{B s , v, q\ sTc \B C , v, 0) = -y/m Bc m Bs Tr (h^I [f^ + a 2 £ T H^) , (3.4) 

where the matrix T multiplies Hy° b ^ on the c quark index because of the c quark spin 
symmetry, and Oi + ao is the most general possible scalar matrix, and multiplies the 
light quark index. In the limit m& ^> m c ^> Aqcd, the addition of momentum q does not 
change the velocity of the meson. The scale of variation of the form factors is controlled by 
the size of the B c bound state. The matrix element Eq. (|3.4| ) measures the overlap of the 
c quark distribution in the B c , the s quark distribution in the B s , and e iq ' x . Equivalently, 
it measures the overlap of the s quark distribution in the B s with the c quark distribution 
in the B c shifted by momentum q. The width of the momentum distribution of the s 
quark is Aqcd, and the width of the momentum distribution of the c quark is of order the 
inverse radius of the B c bound state, 1 ^> Aqcd- Thus a shift in the c quark momentum 
distribution by an amount gCa^ 1 does not affect the overlap amplitude. Consequently, 
the scale of variation of the form factors is Uq , not Aqcd- 

2 The Cl± form factor is the electric coupling, and the ^2 form factor is the magnetic coupling. 
The f^2 form factor gives a B* — > Bj transition amplitude that is not suppressed by powers of 
l/ms, and corresponds to a light quark magnetic moment transition in a quark model. 
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4. The Zero Recoil Limit 

The weak currents sTc, dTc, bTu and bTc do not generate symmetries of the effective 
theory, so their matrix elements at zero recoil cannot be normalized by symmetry consid- 
erations. However, in the limit 1/ao 3> Aqcd i n which the B c is pointlike on the hadronic 
scale Aqcd, it is possible to calculate these matrix elements at zero recoil in terms of 
heavy-heavy bound state wavefunctions and the meson decay constants fs and fo- 

We begin by considering semileptonic B c — > B s decay. In the — > oo limit, the 
kinematics of B c — > B s decay is analogous to that of neutron /3-decay in that the entire 
energy of the decay is transferred to the lepton system. In the rest frame of the B c , 
k° = Ei + E v = m(B c ) — m(B s ), where k is defined in Eq. (|2.1| ), or equivalently q° = 0. 
(The recoil energy of the B s is of order g 2 /m&.) Thus the hadronic form factors only 
depend on g, and the zero recoil point is q = 0. 

The calculation of the form factor at zero recoil proceeds as follows. The initial B c 
state is written asi 



\B cl v) m = J d 3 x $(x) 



(*)^*75^6(-)(0) 



|0> , (4.1) 



where bi \o) creates a b quark with velocity v at the origin, ci + \x) creates a c quark with 
velocity v at the point x, and \l/(x) is the wavefunction of the B c . The superscripts (+) 
and (— ) refer to the f eigenvalue, 

M +) = +ci + \ M _) = -<&~\ (4-2) 



and similarly for bi^ . To compute the Oi form factor in Eq. for B c — * B s semileptonic 



decay, we consider the matrix element of the vector current between the B c and B s at finite 
three momentum transfer q. 



M»(q) = J d 3 z e^ m (B s ,v\V»\B c ,v) nQ 

= J d 3 z e t<r ^ m (B s ,v\s(z)^c v (z)\B c ,v) m 



(4.3) 



where q is the three-momentum transfer to the leptonic system in the decay. Inserting 
Eq. (|4.1|) into Eq. ([Q|), and using heavy field contractions, 

(0\ Cv (x)ci+\y)\0) = ^8(x-y), (4.4) 



3 States with a subscript HQ are normalized to v° rather than to IE. 
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yields 

M»(q) = ijd 3 x e^Vix) n Q (B s ,v\s(x)^ l5 b v (0)\0) . (4.5) 
Performing an operator product expansion 

s(x)^ 5 b v (0) = s(0) 7 ^ 7 5&„(0) + x fc <9 fc s(0) 7 ^75&„(0) + ... , (4.6) 

and using the definition 



(0\s^ l5 b\B s ,v) m = ifB s m B y/y/2m Bs1 (4.7) 
of the B a meson decay constant gives 

^(aoq) = -±=f Bs V^ J d 3 x j* s V(x) , (4.8) 



where we have retained only the first term in the operator product expansion, Eq. 
The above computation also applies to the case where the final meson is Ba instead of 
B s , with the replacement /s sV /^i7 fB d ^/™B~ d - In the SU(3) V limit, fB d ^/™B~ d = 
/b 3 y/rfiB^- The leading ,Si7(3)y-violating correction to this result can be found in ref. 0. 

In the rrib ^> m c ^> Aqcd limit, the wavefunction ^>(x) for the B c is a Coulomb 
wavefunction, so that the form factor Q± can be computed explicitly, 

„ , ^ 1 , , STi 1 / 2 aJ 2 , t s 

^(a q) = —f Bm ^—g—. (4.9) 

There will be corrections to the Coulomb form of the B c wavefunction, because confinement 
effects are significant at the c quark mass. Confinement effects in the meson wavefunction 
at the c quark scale have been studied in detail in the ip system. A quark model with a 
modified Coulomb potential provides a very good description of the spectrum and radiative 
decays for the ip. A similar computation for the B c should provide a good description of 
the B c wavefunction for use in Eq. ( [4.8| ). 

The c and b quark fields were treated as free fields in the computation of Eq. Q4.8Q . 
Radiative gluon corrections are in principle important. Gluon exchange between the b 
and c quarks of the B c and between the b and s quarks of the B s has already been 
included exactly in the definition of the states. The only gluon contributions that are not 
included are radiative gluon corrections where gluons are exchanged between the B c and 
B s . This gluon exchange leads to a violation of factorization in the computation of the 
decay amplitude. The B c does not couple to gluons in the limit that its radius becomes 
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zero. For a finite radius ao, the leading gluon coupling of the B c is to a two gluon operator 
with coefficient proportional to a(j. This produces a small (and incalculable) correction to 
the decay form factors. 

The decay B c —> D proceeds through the quark decay b —>■ u. The D meson is light 
compared with the energy released in the decay, and so can have a large recoil momentum. 
The approximation methods used in this paper cannot be applied in this case. However, 
there is a region of phase space near zero recoil where q < m c , where the heavy c quark 
expansion is still valid. The computation of the Ei form factor in this region is almost 
identical to the computation described above for Q±, with the result 



£i(a <f) = -^=f D V™D I d 3 x e*«-*V(x) 



(4.10) 



Eqs. Q4.8| ) and (|4.10|) both depend on the wavefunction of the B c meson. One can obtain a 
more reliable extraction of KM mixing angles by considering the ratio £i(0)/fii(0), which 
should be insensitive to the detailed form of the B c wavefunction, and thus provides a way 
of measuring and 

The decay B c — > r\ c proceeds through the quark decay b — > c. As for B — > D, the 
heavy quark expansion is only valid in the region q < m c near the zero recoil point. The 
form factor A is calculable in terms of the wavefunctions of the be and cc bound states. A 
straight-forward derivation yields 



A(a Q q) = 2 J d 3 xe-^ 2 % c (x)V Bc (x), 



(4.11) 



where the convolution of the two wavefunctions depends on the radii ao and a v of the B c 



and i] c . The details of this computation are nearly identical to those found in ref. [|10 



for the semileptonic decay of baryons containing two heavy quarks. In the limit that both 
states are described by a Coulomb wavefunction, 



A(a Q q 



3/2 3/2 

Q o a n 
(a + a v y 



1 + 



-»2 2 2 

4(a + a v ) 



(4.12) 



5. Corrections 

Corrections to the results of the previous sections can be divided into two categories. 
The first set of corrections are corrections to heavy quark spin symmetry, and affect the 



10 



relations derived in Sect. 2. The second set of corrections relate to the validity of factor- 
ization and the operator product expansion used in Sect. 4, as well as to the details of the 
B c wavefunction. 

Let us first consider corrections to heavy quark spin symmetry. In this paper, we have 
worked in the limit m& ^> m c ^> Aqcd- There are corrections to the heavy quark theory 
due to the finite mass of the b quark, which are of order Aqcd/^& and m c /mb- These 
corrections are small and will not be discussed further. In addition, there are violations 
of the c quark spin symmetry in the B c . These arise from interactions of the c quark spin 
with the b quark spin, with the orbital angular momentum of the c quark, and with light 
degrees of freedom. The c-b spin-spin interaction is a 1/mj, effect, and is small. There is 
no c quark spin-orbit interaction for the B c because the c quark is in an s-wave. At lowest 
order, the B c is made up of a b quark and c quark in a bound state. There are, however, 
corrections to this form in which the B c wavefunction also contains additional light degrees 
of freedom. In a bag model, this would correspond to exciting gluonic excitations in the 
bag. The spin coupling of the c quark to these light degrees of freedom violates the c 
quark spin symmetry. The interaction energy is of order AQ CD /m c . Since the energy cost 
of exciting a light degree of freedom is of order Aqcd, the net spin symmetry violation in 
the matrix element is of order Aqcd/^c- 

The results of Sect. 4 depend not only on taking m c 3> Aqcd, but also on the size 
of the B c being much smaller than Aq CD . The higher derivative terms in the operator 
product expansion, Eq. ( [Of ) produce corrections of order aoAQCD, because each factor of 
d on the light quark operator produces a factor of Aqcd in the matrix element, and each 
factor of x produces a factor of the size of the bound state a®. There are also violations 
of factorization in the operator matrix element Eq. (|4.3| ). As discussed in the previous 
section, gluon interactions with the B c are of order a;j. There is no suppression factor 
for the interaction of gluons with the B s , since the B s has a size of order Aqcd- Thus 
the gluon interactions produce corrections of order (AqcdOo) 3 - For the B c — > rj c ,xjj decay, 
there is an additional suppression factor of a 3 for the interaction of gluons with the rj c ,ijj, 
so that the net correction is of order Ag^a^a 3 . There are non-perturbative corrections 
to the Coulomb wavefunction of the S c , which are of order Aqcd/^c- As discussed in the 
previous section, most of these effects can be included by modeling the B c by a realistic 
potential which is adjusted to correctly reproduce the B c excitation spectrum. There are 
also radiative corrections which produce corrections of order ct s (mb) and a s (m c ). Finally, 
there are l/m c recoil corrections for B c —>■ D and B c —>■ rj c , ip which are of order q/m c . 
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6. Conclusions 



The semileptonic decay B c — > Diu provides a way of extracting the weak mixing angle 
\V u b\- Theoretical uncertainties can be minimized by extracting the ratio |K i &/V r cs | using 
the ratio of the B c — > D and B c — > B s form factors near zero recoil. The number of 
BcS produced in hadron collisions is much smaller than the number of S's. Nevertheless, 
the B c meson still provides an alternative measurement of \V u b\ to the value which will 
be obtained by comparing semileptonic B — > plv and D — > plv decays. Both of these 
extractions will have corrections due to the finite mass of the c quark whose numerical 
importance will have to be determined experimentally. 
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